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Abstract 

As is well known, hard-sphere crystals of the fee and hep type differ very little in their thermo- 
dynamic properties. Nonetheless, recent computer simulations by Pronk and Frenkel indicate that 
the elastic response to mechanical deformation of the two types of crystal should be quite different 
Q]. By invoking a geometrical argument put forward by R. Martin some time ago Q|, we suggest 
that this is largely due to the different symmetries of the fee and hep crystal structures. Indeed, 
we find that elastic constants obtained by means of computer simulations for the fee hard-sphere 
crystal can be mapped onto the equivalent ones of the hep crystal to very high accuracy. The same 
procedure applied to density functional theoretical predictions for the elastic properties of the fee 
hard-sphere crystal also produces remarkably accurate predictions for those of the hep hard-sphere 
crystal. 

PACS numbers: 62.20.Dc, 82.70.Dd 
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In a recent publication [lj], Pronk and Frenkel reported on a computer simulation study of 
the elastic properties of fee and hep crystals of hard spheres. They found the various elastic 
constants to differ by up to twenty per cent, despite that the thermodynamic properties 
of both types of hard-sphere crystal are almost indistinguishable. Indeed, the free energy, 
pressure and compressibility of the two crystal types deviate from each other by less than 
0.1%, at least for conditions not too far from the melting point In this Comment, we 
point out that the difference in the elasticity of fee and hep hard-sphere crystals is less 
surprising than claimed by Pronk and Frenkel, and that it can be explained by the geometry 
of the packing of the particles within each lattice. 

The relation between the elastic moduli of hep and fee crystals have been studied theo- 
retically and experimentally by number of authors. See, e . Of particular 
interest is the work of R. Martin 4], who derived an approximate transformation of the 
elastic moduli of the fee crystal to those of the hep lattice, making use of the fact that 
both lattice types can be constructed from tetrahedral units. The tetrahedral blocks in the 
fee lattice are equally oriented, while the hep lattice can be built up from two tetrahedra, 
oriented differently from each other and from those of the fee lattice. The transformation of 
any tensor in the fee system of coordinates (defined as usual along the cubic axes) to either 
of the two representations of this tensor in the trigonal geometry of the hep crystal can be 
made by two simple rotations and R^ 2 \ where 
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This suggests that the transformation of the elastic moduli tensor in the fee geometry, 



C FCC , to that of the trigonal geometry of the hep lattice, C HCP , 
average (superposition) of the two trigonal tensors p, 



could simply be the 
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where the subscripts have their usual meaning. It so happens, however, that the two unequal 
tetrahedra of the hep lattice are not independent, but attached to each other. Hence, the 
elastic response of the hep lattice to an external strain should be the combined response 
of both tetrahedra, not just a simple average, implying that Eq. (J2J) requires a correction 
for the internal strain that the connectedness of and interaction between the tetrahedra 
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produce. We refer to ref. 0| for further details. The resulting strain-corrected expression 
for the elastic moduli tensor of the hep lattice, C HCP , reads |4( 



r HCP_ r HCP , HCPyl 

^ijkl — ^ijkl <-^ijrs \y rs tu ) ^tukl , W) 



where C is given by Eq. (j2J) and A is a correction tensor identical to it, except that the 
difference between the two components is taken instead of the sum. 

There are six distinct elastic moduli in the trigonal representation of which only three 
are independent in the fee crystal and five in the hep crystal. The relations between the 
elastic moduli of the fee lattice and those of the hep lattice not corrected for the internal 
strain are given by (using standard Voigt notation) 

CT = (C[r + C[ 2 CC + 2Cl GC )/2 , 
CT = + 5C[ 2 CC - 2C?f c )/6 , 

CT = (C[r + 2C™ G - 2C™ G )/3 , 

CT = {C™° - C[ 2 CG - 2C(n/^ , (4) 
CT = (C[r + 2C™ C + AC™ c )/3 , 
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The strain-corrected constants obey 
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c^ F = cT + (cf 4 OP ) 2 /[jPfr - cT)] . 

The mapping of C FCC onto C HCP implicit in Eqs. (4) and (5) agrees well with experimental 
data on ZnS, a compound that can crystallize both in an fee and in an hep lattice 4]. In 
fact, the mapping works very well for fee and hep crystals of hard spheres too, as we show 
next. 

In Figure 1, we have plotted the relative difference between the various moduli of the fee 
and hep crystals of hard spheres as a function of the dimensionless crystal density ps<J 3 , with 



ps the number density and a the diameter of the spheres, using the computer simulation data 
of Pronk and Frenkel (obtained from Table I of Ref. [|) and the prediction of Martin, Eqs. 
(HHi3) 4]. We find that the largest deviation between the two is about 6%. This implies that 
the approximate theory outlined above, based entirely on a geometric argument, plausibly 
explains the difference between the elastic moduli of fee and hep crystals. That is, geometry 
explains almost all of the found differences in elastic behavior of fee and hep crystals. 

In order to further verify Eqs. (JUEJ), we calculated the elastic moduli of the hep hard- 
sphere crystal, using results for the elastic moduli of fee crystals of hard spheres obtained 
from density functional theory (DFT) ~\. We used the predictions of the modified weighted- 
density approximation DFT, MWDA DFT, because they are known to agree very well with 
the results of computer simulations. The results of the mapping of the fee moduli onto the 
hep moduli are presented in Figure 2, again as a function of the dimensionless density of the 
spheres. The agreement with the results of the computer simulations of Pronk and Frenkel 
is quite good for all moduli except C13, for which it is not as impressive but still satisfactory. 

We believe to have demonstrated that the large difference in elasticity between the fee 
and hep crystals of hard spheres found in jjj is largely caused by the geometrical differences 
of these two types of crystal lattice. We point out that even when the thermodynamic 
properties of fee and hep crystals are similar to the point of being virtually indistinguishable, 
there is in fact no reason for their elastic properties to be similar too. The reason is that 
similarities in the free energy landscape at long wavelengths do not preclude differences at 
short wavelengths. Indeed, (direct) correlations at short wavelengths contribute much more 
significantly to the elastic constants than to the equilibrium free energy ^. 
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FIG. 1: Relative difference = \C^ CC - Cf° p \jCff c between the elastic moduli of the hep and 
fee crystals, as a function of the dimensionless density ps«r 3 . Shown are the results of computer 
simulations of Ref. l], 6ij, and the ones computed using relations Eqs. (^SJ), S[j . 
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FIG. 2: The dimensionless elastic moduli Cy = Cija 3 /ksT of the hep crystal computed using the 
relations Eqs. and the results of the MWDA DFT calculations for the fee crystal L 7], as a 

function of the dimensionless crystal density /O5O" 3 . Also shown are the results of the computer 
simulations of Ref. Illl . 
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